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SUMMARY 

The  purpose  of  this  experiment  was  to  measure  the  scaling  with 
distance  of  the  transverse  mode  stimulation  in  the  free  electron  laser 
experiment  at  LURE  in  Orsav.  These  experiments  were  performed*  and  the 
analysis  has  revealed  that  the  phenomenon  scales  as  predicted  by  the 
small  sianal  theory.  A  new  effect  was  observed  with  a  multiple  mode 
incut  laser  probe?  the  theory  t.o  describe  this  effect  was  developed.  In 
further  analysis  of  data  taken  previously*  we  were  also  able  to 
demonstrate  the  scaling  of  the  transverse  coupling  terms  in  the  sain  as 
a  function  of  the  electron  beam  size  and  the  resonance  parameter. 

Our  main  conclusion  from  this  work  is  that  the  size  of  the 
off-diagonal  terms  in  the  transverse  sain  matrix  is  of  the  order  of 
unity  if  the  electron  beam  and  the  laser  beam  sizes  are  comparable. 

Since  this  situation  is  likely  to  hold  for  all  FEL  devices*  transverse 
effects  cannot  in  general  be  neglected.  To  date*  these  effects  have  not 
shown  up  due  to  the  low  gain  of  the  oscillators  which  have  been 
orerated.  For  devices  with  a  aain  of  the  order  of  unity  or  higher. - 
djffractive  effects  will  stron3ly  modify  the  laser  mode.  Programs 
developing  FELs  in  the  high  sain  regime  should  be  aware  of  this  fact* 
and  budget  the  research  time  to  investigate  the  implications  to  their 
particular  configuration. 

This  experiment  is  the  first  capable  of  resolving  an  excitation  of 
t ho  higher  order  transverse  modes  in  a  free  electron  laser*  and  has 
established  the  existence  of  off-diasonal  terms  in  the  transverse  gain 
matrix.  The  technique  we  have  developed  is  the  only  diagnostic  of  the 
transverse  mode  mixing  effect  which  can  diagnose  its  presence  in  devices 


with  optical  gain  below  approximately  10 Q‘i  per  pass 


INTRODUCTION 


The  purpose  oF  this  experiment  was  to  obtain  some  additional 
information  on  the  phenomenon  of  mode  mixing  in  the  FEL recently 
discovered  at  Orsav  [Appendix  13.  u  sinale  transverse  mode  laser 
beam  is  observed  to  produce  stimulated  emission  into  a  series  of 
hiaher  order  modes.  Me  have  developed  a  technique  which  enables  us 
to  measure  this  effect.  In  this  work  ?  we  have  obtained  a  new  set 
of  data  and  have  been  able  to  show  both  qualitative  and  suanti tave 
asreement  with  the  theory  [Appendix  II 3. 

Transverse  mode  mixina  is  important  to  the  operation  of  FEL 
hiah  gain  oscillators  and  amplifiers.  For  low  sain  systems,  the 
vector  of  mode  amplitudes  and  phases  is  little  chansed  in  one 
passi  and  the  power?  which  is  quadratic  in  the  amplitudesr  only 
increases  noticeably  in  the  modes  populated  by  the  input  beam. 
Oscillators  operated  in  this  resime  will  show  no  effects  from  mode 
mixina.  At  moderate  values  of  aain.-  the  rotation  of  the  input  mode 
vector  can  be  substantial?  resulting  in  multiple  mode  output 
regardless  of  the  input  mode.  Storage  ring  lasers  are  expected  to 
orerate  in  this  reaime.  If  the  sain  is  sufficiently  large?  the 
propaaatins  li3ht  wave  can  be  very  strongly  modified.  In  this 
regime?  a  growing  field  distribution  can  be  maintained 
which  is  substantially  better  Focussed  than  the  lowest  order 

Gaussian  mode  of  free  space.  This  situation  is  likely  to  be 
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present  in  the  amplifier.  r 
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Tho  presence  of  cross  terms  in  the  transverse  aain  matrix  was 
First  observed  oxper imental 1 y  a  year  aao  [Appendix  13.  A  basic 
theory  was  derived/  and  applied  to  the  conditions  oF  the 
measurement  system.  This  theory  revealed  scalins  relations 
which  had  not  been  investiaated  in  the  oriainal  experiments.  In  this 
work/  the  experiments  were  extended  with  additional  equipment  to 
enable  us  to  probe  the  scalins  as  a  Function  oF  the  analxzins 
aperture  displacement/  and  the  older  data  was  investiaated  to 
check  the  scalina  with  the  resonance  parameter  and  the  electron 
beam  size.  The  theory  was  also  extended  anal  ytical  1  ■/  to  cover  the 
case  oF  an  arbitrary  combination  oF  modes  in  the  input  beam.  The 
results  show  clearly  that  the  seneral  theory  is  correct  in  the 
small  sisnal  resime.  In  the  larse  sianal  reaion/  where  the 
application  oF  tho  theory  is  still  incomplete.-  the  experimental 
method  we  have  developed  may  prove  useful  as  a  tool  to  identify 


the  effects  of  the  mode  mixina 


DESCRIPTION  OF  THE  EXPERIMENT 


The  transverse  mode  content  of  the  optical  Field  propasatins 
throush  tho  FEL  can  be  conveniently  thousht  of  as  a  vector  whose 
elements  are  the  amplitudes  and  phases  of  the  modes.  As  »he  beam 
ptopaaates  throush  free  sparer  the  phases  evolve  in  a  fashion 
characteristic  of  each  moder  but  the  amplitudes  remain  constant 
(in  the  absence  of  losses).  The  amplifier  contributes  an 
additional  vector  of  complex  field  amplitudes  which  describes  the 
stimulated  radiation/  and  which  is  in  general  not  a  simple 
multiple  of  the  input  vector.  The  Field  vector  at  the  output  of 
tho  FEL  amplifier  is  the  sum  of  the  input  vector  propasated  to  the 
output  and  tho  stimulated  emission  vector.  The  magnitude  squared 
of  the  output  vector  will  be  larger  than  that  of  the  input  by  a 
factor  which  is  the  net  amplification/  and  the  vector  will  be 
rotated  in  the  mul t i -d imens i onal  mode  space.  If  the  sain  is  larse 
and  mode  mixing  occurs/  the  output  field  vector  will  be  strongly 
rotated.  Under  these  o i rcums tances r  the  laser  can  be  expected  to 
run  multi-mode-  and  a  self-consistent  calculation  of  the  field 
distribution  must  be  performed  to  identify  the  steady  state  field 
vector  and  the  gain  it  experience?.  In  a  low  sain  oscillator/ 
however/  the  same  level  of  mode  mixing  will  produce  no  effects 
because  the  relation  of  the  vector  is  negligible. 

The  gain  of  the  Orsax  FEL  is  very  low/  on  the  order  of  10A-3. 
Even  though  the  conditions  for  mode  mixing  are  present  (the 
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dimensions  of  laser  and  electron  beam  are  comparable) >  the 
rotation  of  the  field  vector  is  minuscule.  As  could  be  expected 
under  these  conditioner  the  oscillator  produced  a  near) y  pure 
TEMoo  Gaussian  mode  Cl  3.  < Even  the  Stanford  oscillator  which  has  a 
sain  of  about  107£  C23  Produces  a  Gaussian  mode  C33).  The  newer  FEL 
systems  now  under  construction  at  Stanford/  Orsa/.-  and  elsewhere 
are  projected  to  have  sains  ransins  up  to  a  factor  of  ten  per 
pass.  The  mixing  phenomena  of  the  transverse  modes  will  play  an 
important  role  in  the  operation  of  these  new  devices.  Our 
technique  permits  us  to  observe  the  three  dimensional 
amplification  process  even  in  devices  whore  the  sain  is  so  low 
that  no  effects  can  be  observed  in  the  modes  of  the  oscillator. 


By  means  of  an  adjustable  aperture  placed  in  the  amplified 
beam/  we  have  been  able  to  measure  the  excitation  of  higher  order 
transverse  modes  induced  by  the  nonuniform  transverse  profile  of 
the  electron  beam.  This  information  is  mapped  into  the  spatial 
distribution  of  the  electric  field  at  the  analyzing  aperture.  The 
ratio  of  the  transmitted  stimulated  power  to  the  transmitted  power 
of  the  externa]  laser  is  measured  as  a  function  of  the  aperture 
diameter.  The  functional  dependence  of  the  ratio  on  the  aperture 
diameter  is  a  diagnostic  on  the  transverse  mode  mixing  which  has 
occurred  during  the  amplification  process. 


The  apparatus  we  used  is  essentially  identical  to  that 
reported  in  the  original  3ain  measurements  C4-73.  A  chopped 
external  argon  laser  is  Focussed  coaxially  onto  the  electron  beam 
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in  the  optical  Klystron  where  it  is  amplified  by  the  FEL 
interaction.  The  output  beam  is  passed  through  an  adjustable 
aperture*  filtered  in  a  monochrome t er »  and  detected  in  a  resonant 
detector.  The  amplified  power  is  extracted  in  a  double 
demodulation  scheme:  a  13  Mhz  locK-in  detects  the  sum  of  the 
spontaneous  and  the  stimulated  power*  and  a  subsequent  lock-in 
operatins  at  the  frequency  of  the  chopper  extracts  a  signal 
proportional  to  the  stimulated  power  at  the  detector. 

In  this  experiment*  the  undulator  sap  is  scanned  to  set  a 
maximum  of  the  optical  Klystron  sain  spectrum  £71  at  the  laser 
wavelength.  Alignment  is  performed*  and  the  powers  of  the 
transmitted  laser  and  of  the  stimulated  emission  are  measured  in 
separate  channels  and  recorded  for  about  100  seconds.  The  laser  is 
then  blocked  in  order  to  permit  a  measurement  of  the  noise 
background!  the  aperture  diameter  is  readjusted  and  realigned  (if 
necessary).1  and  the  two  power  measurements  are  repeated.  In  the 
absence  of  cross  terms  in  the  gain  matrix*  the  ratio  of  the  two 
ch£innels  is  independent  of  the  aperture  diameter.  Excitation  of 
higher  order  modes  is  signaled  by  a  deviation  of  this  ratio  from 
unity.  The  detection  system  was  checked  for  spurious  sources  which 
might  mimic  this  signal;  any  deviation  from  linearity  or  position 
sensitivity  of  the  detector  was  verified  to  be  small  compared  to 
the  observed  signal. 

The  interpretation  of  the  results  demands  an  accurate 
Knowledge  of  the  eleotron  and  optical  beam  sizes.  Those  were 


measured  with  the  aid  of  a  Reticon  diode  array  of  25  micron 
spac ina»  and  the  digitizing  oscilloscope.  The  shape  and  the  full 
width  at  half  maximum  of  the  input  laser  mode  was  recorded  at  five 
positions  about  the  Focus.  both  in  the  vertical  and  the  horizontal 
dimensions.  Unf ortunatel y .  the  Arson  laser  provided  by  Stanford 
For  this  worK  is  now  aseins  noticeably,  and  it  was  impossible  to 
obtain  sinsle  transverse  mode  operation.  In  order  to  see  the 
effects  of  the  multiple  mode  input,  a  seneral  reformulation  of  the 
theory  was  done,  including  an  arbitrary  number  of  modes  in  the 
input  beam  (see  Appendix  II.  equations  *  1  —  1 0  >  > .  From  this 
analysis,  it  follows  that  an  analytic  comparison  with  the  theory 
would  require  a  Knowledge  of  the  amplitudes  and  the  relative 
phases  of  all  the  modes  present  in  the  input  beam,  an  impossible 
tasK  with  the  available  instrumentation.  On  the  other  hand, 
qualitative  comparison  with  the  theory  still  reveals  the  presence 
of  several  important  trends  in  the  data,  as  we  shall  see  below. 

The  electron  beam  size  is  measured  with  the  aid  of  the 
synchrotron  radiation  emitted  from  the  beam  in  the  dipole  magnet 
in  Front  of  the  optical  Klystron.  The  vertical  and  the  horizontal 
beam  shapes  ar«  typically  observed  to  be  Gaussian,  with  a  width 
which  varies  according  to  the  current  in  the  ring.  These 
dimensions  also  change  with  the  tunes  of  the  ring,  so  separate 
electron  beam  measurements  are  required  on  each  injection  during 
the  experiment. 

Due  to  a  technical  problem  with  the  extraction  mirror,  the 
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experiments  were  performed  this  time  with  the  probe  laser  beam 
exiting  the  vacuum  system  through  the  rear  oscillator  mirror.  The 
transmission  of  this  mirror?  a  hish  reflector  at  6328  A t  is  about 
80%  at  the  Arson  514?  A  wane  1 enath .  It  acts  as  a  defocussins  lens 
with  a  focal  lensth  of  about  -5.2  m.  These  experiments  also  had  to 
be  performed  at  low  onersv  •  <  i=6  HeO)  in  order  to  avoid  damasins 
C81  the  oscillator  mirrors  with  the  UO  senerated  in  the  optical 
Klystron.  This  produced  a  somewhat  lower  s i anal -t o-no i s e  ratio 
than  in  the  previous  worK. 


nP=ER"ATIQM  OF  HIGHER  ORDER  MODE  GENERATION 


As  reported  in  Appendix  II,  we  completed  a  set  of 
measurements  of  the  sain  as  a  function  of  the  iris  aperture  for 
two  distances  of  the  aperture  from  the  optical  Klystron*  one  as 
close  as  possible  <2.9  m )  and  the  other  essentially  at  infinity 
(14.3  ml.  These  measurements  show  the  evolution  of  the  phases  of 
the  individual  components  predicted  in  Appendix  I.  Measured  near 
the  Klystron  'Fisure  7  o'-'  Appendix  1 1  >  r  the  phases  have  shifted 
little  from  their  values  in  the  interaction  reaion.  Measured  at 
infinity  (Fisure  <?  of  Appendix  III.-  the  phases  have  attained  their 
maximum  valuer  depressins  the  sain  measured  at  the  center  of  the 
beam  relative  to  that  of  the  near  :one  measurement.  The  results 
exhibit  the  qualitative  behavior  expected  from  the  theory. 

Our  data  analysis  proaram  also  addressed  the  data  from  the 
earlier  experiments  of  a  year  aso,  Quantitative  asreement  has  now 
been  shown  b.tween  the  predictions  of  the  theory  and  the  results 
of  the  earlier  sinslo  mode  experiments.  An  examination  of  f isure  4 
of  Appendix  II  confirms  beyond  a  doubt  that  the  low  field  theory 
is  understood. 

The  theory  also  predicts  that  the  sain  matrix  should  be 
essentially  independent  of  the  resonance  parameter.  We  were  able 
to  taKe  advantase  of  the  optical  Klystron's  unique  spectral  shape 
to  verify  »ho  scalina  of  the  mode  mix  ins  as  a  function  of  the 
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manetic  oap.  Figure  6  of  Appendix  II  fhous  three  scans  taken  at 
different  ‘.'a lues  of  the  resonance  parameter  and  superimposed. 

These  three  scans  correspond  to  gain  and  absorption  peaks  four 
klystron  f  ri  rises  away  from  resonance.  They  agree  within  the 
experimental  errors. 

Figures  5  and  8  demonstrate  the  scalina  of  the  effect  as  a 
function  of  the  electron  beam  size.  These  f  isures.-  taken  with  both 
single  and  multiple  mode  input',  show  unmi s t aKeab 1 y  that  laraer 
electron  beams  show  less  excitation  of  the  higher  order  transverse 
modes.  This  trend  is  expected  because  in  the  limit  of  a  large 
electron  beam,  the  gain  medium  exactly  reproduces  the  input  mode 
combination  so  that  no  aperture  dependence  will  be  observed. 

This  array  of  data  completes,  in  our  view,  the  verification 
of  the  small  signal  theory  of  transverse  mode  stimulation.  In 
future  systems  which  operate  with  high  gain  or  in  the  saturated 
regime,  one  can  anticipate  that  novel  effects  (such  as  "sain 
Focussing")  will  appear.  We  can  now  approach  those  systems  with 
the  assurance  that  the  Fundamental  equations  are  well  known 


and  verified. 
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Abstract.  We  derive  the  most  general  equations  of  motion  for  the  electrons  and  the 
electromagnetic  field  in  a  free-electron  laser  including  the  effects  of  diffraction  and  pulse 
propagation.  The  field  evolution  is  expressed  in  terms  of  the  amplitudes  and  phases  of  a 
complete  set  of  transverse  modes.  The  analytic  solution  is  given  in  the  small-signal  regime, 
where  the  theory  is  shown  to  be  in  excellent  agreement  with  a  recent  experiment  at  Orsay. 

PACS:  42.60,  42.20,  42.55 


Stimulated  by  the  original  free-electron-laser  experi-  the  minimum  number  of  physically  observable  quan- 

ments  [1]  in  1977,  a  number  of  authors  have  contri-  tities:  the  transverse  optical  modes  of  the  system.  The 

buted  to  the  development  of  a  purely  classical  theory  field  evolution  is  expressed  in  terms  of  a  complete  set  of 

for  the  electron  dynamics  and  the  electromagnetic  orthogonal  transverse  modes;  equations  are  developed 

wave  growth  in  these  devices.  The  initial  work  as-  for  the  propagation  of  the  amplitude  and  phase  of  each 

sumed  the  light  could  be  represented  by  a  single-  mode.  In  physical  systems  which  operate  on  a  few  of 

frequency  plane  wave  [2,  3].  The  first  generalization  the  lowest-order  modes,  this  approach  greatly  in- 

was  required  to  explain  the  extremely  short  pulse  creases  the  accuracy,  and  may  reduce  the  required 

phenomena  observed  at  Stanford  [4,  5].  The  inclusion  computer  time  for  the  calculation  by  working  in  a 

of  the  longitudinal  modes  in  the  theory  [6-8]  per-  vector  space  well  matched  to  the  solution  of  the 

mitted  the  explanation  of  the  cavity  detuning  curve,  problem.  For  the  oscillator  case,  the  appropriate 

and  predicted  a  range  of  phenomena  in  the  pulse  choice  of  modes  is  the  set  of  eigenmodes  of  the  cavity, 

structure  which  have  yet  to  be  observed.  More  re-  For  the  amplifier,  the  vector  space  of  modes  is  de¬ 
cently,  the  theory  has  been  broadened  to  include  the  termined  by  the  characteristics  of  the  input  mode, 

transverse-mode  structure  of  the  optical  beam  [9-15].  which  is  presumably  a  TEMno  Gaussian  mode.  In 

Until  our  work  at  Orsay  [16],  no  experimental  infor-  either  device,  an  optimum  design  would  result  in  the 

mation  has  been  available  to  test  the  validity  of  these  excitation  of  as  few  of  the  higher-order  modes  as 

so-called  3D  theories.  possible.  The  modal  decomposition  method  is  there¬ 
in  this  paper,  we  present  a  new  approach  to  calculating  fore  well  adapted  to  the  prediction  and  optimization  of 

the  three-dimensional  effects  operative  in  free-electron  the  operation  of  the  free-electron  laser  (FEL). 

lasers.  The  previously  mentioned  approaches  consider  In  the  first  section,  we  derive,  in  their  most  general 

the  growth  of  the  field  6(r,  r)  along  the  propagation  or  form,  the  equations  governing  the  dynamics  of  the 

z  axis  by  evaluating  its  change  at  each  point  (.x,y),  and  complex  mode  amplitudes.  The  subsequent  sections 

integrating  numerically  through  the  interaction  region  reduce  these  equations  to  the  familiar  case  of  the  small 

in  the  time  domain  [9-14]  or  in  the  frequency  domain  signal,  low-gain  result  (Sect.  2).  Here,  the  problem 

[15].  These  techniques  all  demand  long  computer  runs  becomes  linear,  the  mode  evolution  can  be  described 

if  they  are  to  be  applied  to  a  real  experimental  by  a  matrix  transformation,  and  we  retrieve  the  well 

situation.  Our  approach  decomposes  the  problem  into  known  gain  equation  complete  with  filling  factor.  The 
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theory  is  then  applied  to  the  case  of  the  Orsay  experi¬ 
ment.  where  the  results  are  in  excellent  agreement  with 
an  experiment  [16]  performed  recently  which  exhibits 
the  off-diagonal  terms  of  the  gain  matrix. 


I.  Theoretical  Development 
of  the  Fundamental  Equations 

The  FEL  system  is  properly  described  by  the  coupled 
Maxwell  and  Lorentz  force  equations.  From  these,  we 
shall  derive  a  self-consistent  set  of  equations  describing 
the  electron  and  the  transverse  optical  mode  dynamics. 
We  use  the  dimensionless  notation  originally  deve¬ 
loped  by  Colson  (in  fact  this  work  is  a  generalization  of 
Colson’s  work  to  include  transverse  modes  and  we 
shall  stay  as  close  as  possible  to  his  original  notation). 
Let  us  recall  his  main  equations  describing  the  field 
and  electron  dynamics  in  the  slowly  varying  phase  and 
amplitude  approximation  [18]: 

^  =  a'  cos(C  +  <f>),  (1) 


dx  V' 
da 

—  =  -,'<e-%.vo, 


where 

C{t)  =  (k  +  k0):(t)-cjt 

is  the  dimensionless  electron  phase, 

v(t)  =  L[(k  +  k0)  P.(t)  -  fc] 

the  dimensionless  resonance  parameter. 


Here  we  consider  an  /V-period  helical  undulator  of 
length  L,  magnetic  period  /0  =  2rr  peak  magnetic 
field  B,  and  deflection  parameter  K  =  93.4B 
[Gauss]  A0  [cm].  An  electron  beam  of  energy  /me2, 
and  number  density  q  travels  along  the  axis  of  the 
undulator:  an  individual  electron  has  longitudinal 
coordinate  r(t)  and  longitudinal  velocity  c/L(r)  at  time 


t.  A  helically  polarized  plant’  wart’  of  wavelength 
/.  =  2r.  k.  frequency  o.  and  electric  field  fiiz.n  -  Elz.t) 
exp{i[kr-wr  +  d>(r. r)]J  interacts  with  the  electrons.  In 
(3),  <  is  the  average  over  the  initial  phase  and 
resonance  parameter  v0  of  the  electron  population  at 
the  position  z. 

Equations  (1,2)  are  derived  directly  from  the  Lorentz 
force  equation  and  describe  the  effect  of  the  radiation 
field  on  the  electrons.  The  work  done  by  the  longitu¬ 
dinal  field  on  the  electrons  is  neglected  here,  which  is  a 
good  approximation  provided  that  the  modes  are  not 
too  divergent  [14]  /./oj0<£2n2KN/y.  Equation  (3)  is 
derived  from  the  Maxwell  equations  and  describes  the 
effect  of  the  electron  on  the  radiation  field.  The  set  (1), 
(2),  and  (3)  is  self-consistent.  Indeed,  those  equations 
are  very  close  to  being  the  most  general  classical 
equations  describing  the  FEL  dynamics.  They  apply  to 
high  and  low  gain  devices  (r  >  1  or  r<?  1).  high  field 
and  low  field  cases  (u'  £>  1  or  a'  <?  1 ),  and  include  the 
effects  of  multiple  longitudinal  modes  (laser  lethargy 
effects)  through  the  z  dependence  of  r ,  E,  and  <£.  Slight 
modifications  allow  their  extension  to  the  cases  of: 

-  the  planar  undulator  [19], 

-  the  tapered  undulator  [18], 

-  the  optical  klystron  [14],  and 

-  space  charge  effects  [19]. 

However,  the  plane-wave  approximation  cannot  ac¬ 
curately  describe  the  transverse  effects  produced  by  the 
finite  transverse  extent  of  the  optical  mode  and  the 
electron  beam.  A  filling  factor  calculated  with  an 
ad-hoc  overlap  integral  can  be  added  to  the  results  of 
this  calculation,  and  gives  satisfactory  results  in  the 
small  signal  regime  only  so  long  as  one  is  not  interested 
in  the  exact  transverse  field  profile. 

To  relieve  this  last  restriction  on  the  theory,  we  assume 
the  field  to  be  described  in  free  space  by  the  paraxial 
wave  equation  [20] : 


£(r,  f)ei*lr'"  =  0. 


(9) 


This  equation  is  derived  from  the  wave  equation 
(V2-c~2c2/tV)<?  =0  in  the  slowly  varying  amplitude 
and  phase  approximation,  and  has  been  widely  used  in 
laser  field  calculations  [  1 7.  20].  The  general  solution  of 
(9)  can  be  expressed  as  a  linear  combination  of  a 
complete  set  of  orthogonal  modes.  If  we  define  these 
modes  by  the  complex  amplitude  E„c\ p(ivm).  where 
Em  is  real  and  in  is  the  generalized  index  of  the  mode  (in 
the  two-dimensional  transverse  space  we  consider,  in 
represents  two  integer  numbers),  the  most  general 
expression  for  the  field  is 


E(r.fle,*""  =  YcJt)Ejr)c"-''\  (10) 

where  cm  is  complex  and  time-independent  in  free 
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space.  The  orthogonality  relation  reads 
”£e'v-  =  ()  . 


f^E.e- 


jnv; 


(11) 


where  we  have  chosen  a  convenient  normalization 
which  makes  the  Em  dimensionless.  The  modes  can  be 
chosen  in  a  variety  of  symmetries,  but  it  is  useful  to 
exhibit  their  specific  form  in  cylindrical  symmetry : 


£„i(r)  = 


<Mr)= 


1/  2'  *  2p  i  *0  /  r  V  fee 

!  (p  +  /) !  ( 1  +  <50l)  >v(r)  \w(r)/  \si 


cos/01 
sin/OJ 


(12) 


kr 2 


<13) 

where  r  is  the  radial  and  9  is  the  azimuthal  coordinate, 
w0  is  the  beam  waist,  L‘p( 2r2/w2)  is  the  associated 
Laguerre  polynomial,  and 

w2(z)  =  wg(l  +  ^j),  (14) 


,2 

-0 


*(--)  =  r  1  + 


,2 

-0 


Jtw; 


(15) 


(16) 


These  modes  are  very  useful  for  the  case  of  a  cylindri¬ 
cal  electron  beam  aligned  to  the  axis  of  the  light  beam. 
For  an  ellipsoidal  electron  beam  profile,  or  off-axis 
electron  injection,  the  rectangular  eigenmodes  are 
more  appropriate.  Although  we  will  use  the  cylindrical 
modes  in  the  examples,  we  proceed  with  the  general 
theoretical  development  which  makes  no  assumptions 
on  the  specific  form  of  the  modes. 

In  FEL,  the  coefficients  in  (10)  become  time  dependent. 
We  wish  to  calculate  the  evolution  of  the  amplitude 
and  phase  of  these  mode  coefficients.  Proceeding 
through  the  derivation  of  (1-3),  making  only  the  slowly 
varying  amplitude  and  phase  approximation,  but  now 
using  (10)  and  (11),  we  find 

JZ  =  I  laj  cos(C  +  H'„  +  <t>„) , 


vx 

dA. 

dz 

X 


f dxdy  _ 
=  -  J - -  r£_e 


nwz 


”<e  i;>. 

'  '  -.o 


where  we  have  made  the  new  definitions 
4  neSLK 

^  Sn2e2.\L2K2 

>iT,t)  =  — '-3  — - e(r.r). 


WC* 


cjr,f)  =  |cjr,t)|  e,*'"u-°. 


(17) 

(18) 

(19) 

(20) 

(21) 

(22) 


As  before.  (17)  and  (18)  describe  the  effect  of  the 
radiation  field  on  the  electrons,  and  (19)  describes  the 
growth  or  decay  of  the  radiation  field  due  to  its 
interaction  with  the  electrons.  The  change  in  (17)  is 
quite  straightforward.  Equation  (19)  shows  clearly  the 
fact  that  the  growth  in  the  m,h  mode  amplitude  and 
phase  is  given  by  the  operlap  integral  of  the  inphase 
and  out-of-phase  components  of  the  charge  density 
with  the  complex  conjugate  of  that  mode,  as  one  would 
expect.  We  note  that  the  only  assumptions  made  on 
the  modes  £mexp(itpm)  used  in  (17-19)  are  orthogo¬ 
nality  and  completeness.  This  means  these  equations 
are  also  valid  for  the  cases  of  waveguide  modes  and 
dielectrically  loaded  cavities.  In  this  case,  w0  is  no 
longer  the  mode  waist  in  the  usual  Gaussian  sense,  but 
is  defined  by  (11).  As  before,  these  equations  are  self- 
consistent.  An  example  of  this  fact  is  i.  energy 
conservation  equation 


^l\aJ2=-2SdXdy 

dxt  3  Ttw; 


(23) 


which  is  derived  from  (17)  and  (19).  The  left-hand  side 
of  (23),  the  total  energy  gained  by  all  the  modes,  is 
equal  to  the  energy  loss  integrated  over  all  of  the 
electrons  in  the  beam. 

Equations  (17-19)  retain  all  of  the  generality  of  (1-3). 
They  are  valid  for  high  and  low  fields,  and  high  and 
low  gain  systems.  They  take  into  account  the  evolution 
of  the  transverse  modes  explicitly,  and  the  evolution  of 
the  longitudinal  modes  implicitly,  by  keeping  track  of 
the  r  dependence  of  the  charge  density  r(r.  t)  and  of  the 
mode  amplitudes  a(r,  t).  For  simplicity  in  the  follow  ing 
development,  we  drop  the  explicit  i  dependence  which 
has  been  thoroughly  discussed  by  Colson  [18],  and 
concentrate  on  the  transverse  phenomena. 

As  discussed  in  [19],  the  generalization  to  the  case  of 
the  planar  undulator  is  no  more  than  a  change  in  the 
definition  of  the  two  parameters 


(24) 

(25) 


where 


UJ1=J  0 


K2 


4  +  2  K- 


K2 


■ 2I  J'U  +  2K2r 


(26) 


Equations  (17-19)  can  be  integrated  numerically  to 
find  the  evolution  of  the  optical  wave  in  any  Compton 
regime  FEL.  In  a  high-field  experiment.  (18)  and  (19) 
are  nonlinear  in  a.  and  the  wave  evolution  can  only  be 
obtained  numerically.  In  this  ease.  (17-19)  provide  a 
precise  and  efficient  technique  for  solving  the  general 
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problem.  In  a  low-field  situation  such  as  we  find  at 
Orsay,  however,  the  problem  becomes  linear,  and  can 
be  solved  analytically.  We  proceed  with  the  low  field 
case  in  the  next  section. 


2.  The  Low-Field  Solution 

The  low-field  case  is  defined  by  |ap|  <?  1  for  every  mode. 
In  other  words,  the  electrons  do  not  become  over¬ 
bunched.  Experiments  which  operate  in  this  domain 
include  the  low-field  amplifier  experiments,  and  sto¬ 
rage  ring  FEL  oscillators  which  saturate  by  mecha¬ 
nisms  other  than  overbunching.  The  ignition  of  any 
FEL  oscillator  also  occurs  in  this  domain. 

2.1.  The  Gain  Matrix 

Equations  (17-19)  can  be  solved  by  integrating  (17) 
and  (18)  to  lowest  order  in  the  fields  am  and  inserting 
the  result  for  £  into  (19).  If  the  electrons  are  uniformly 
distributed  initially  in  phase,  we  find 


da^)  r 


-j—  =  J  dz'  1  </t"Mm,(T,  z")ajt") ,  (27) 

UT  0  0 

where 

1  ,'dxdv 

MJi?,  T")  =  X  J  ■'~T  }’)  £JX’  >’  ri  £JX-  y-  T") 

2  n  w0 

.  e  -  y.  t)  -  y.  filing  -  i'o(t  -  j  ->g  J 

Equation  (27)  describes  a  linear  evolution  of  the  mode 
amplitudes,  and  upon  integration,  gives  the  relation 

«m(T=l)  =  (/  +  G)m„un(r=0),  (29) 

where  /  is  the  identity  matrix,  and  G,  which  is  generally 
not  Hermitian,  has  elements 

l  T  < 


9mn=  t") 


0  0  0 


+  j*.  JdTjAfjT,.r3) 

0  0  0 

•J  *4  J  J'  dx6Mln (T,.  r6)  + ....  (30) 

0  0  0 

The  higher-order  terms  in  <y„m  are  proportional  to  r: 
and  higher  powers  of  r,  and  are  negligible  in  the  low 
gain  case. 

Evidently  this  matrix  is  of  great  interest  since  multiple 
passes  of  the  electron  beam  will  result  in  multiple 
products  of  this  matrix,  greatly  simplifying  the  calcu¬ 
lation  of  the  modes'  growth.  We  shall  discuss  the 
consequences  for  an  oscillation  experiment  in 
Sect.  2.2. 


Let  us  note  that  this  gain  matrix  is  generally  complex 
and  defines  the  growth  of  the  amplitude  of  the  field. 
Sometimes  people  speak  of  the  gain  in  a  mode  "in"  as 
the  energy  gained  by  this  mode  in  a  pass  through  the 
undulator.  This  gain  is  simply  2  Re  t  Of 

course,  one  must  keep  in  mind  that  energy  is  radiated 
into  other  modes,  and  that  cross  terms  will  mix  a 
multiple  mode  input.  If  the  input  beam  is  truly  mo¬ 
nomode.  the  power  radiated  into  the  nlh  mode  i;  lower 
than  that  into  the  m'h  mode  by  the  ratio 
lsmJ2/2  Re  which  is  small  for  low  gain  (r  <g  1> 
systems.  It  is  only  in  this  case  that  it  makes  sense  to 
speak  of  the  gain  of  a  mode.  In  high-gain  systems, 
however,  the  off-diagonal  terms  can  lead  to  substantial 
emission  of  energy  into  the  higher-order  transverse 
modes.  If  the  input  beam  is  multimode,  of  course, 
mode  mixing  occurs  at  all  power  levels. 

Let  us  calculate  gm„  in  the  simple  case  of  experimental 
interest  where  the  electron  beam  is  cylindrical,  and  a 
good  choice  of  modes  is  the  cylindrical  cavity  eigen- 
modes  (12)  and  (13).  We  restrict  ourselves  to  the 
weakly  diverging  case  nwl$>/.L  where  the  gain  takes 
on  its  most  familiar  form.  The  mode  amplitudes  and 
phases  in  (28)  become  independent  of  r.  and  we  can 
integrate  the  first  term  in  (30)  to  find  the  gain.  The 
average  over  the  resonance  parameter  in  (28)  becomes, 
under  the  assumption  of  a  Gaussian  distribution  of 
centroid  v,  and  deviation  o . 


<e =  e  -  e-,,‘,r~;  l.  (31) 

Under  the  weak-divergence  approximation,  and  as¬ 
suming  negligible  pulse  slippage  effects  (long  electron 
bunch  length  a,  >/V/.),  the  only  time-dependence  in 
(28)  is  that  of  (31).  For  small  spread  av  <S  1  the  integral 
gives  the  well  known  gain  spectrum 


gmn-  j— ; r'ix.y)  EJx.y)  £„|.v,y)  e'iVm<x,>'’  eiv"u,!” 

7TH  n 


1  —  cost;.  —  sin  i\ 


v ,  v 

-  v  cosv.  +  sin  v 

i  i  * 


In  the  usual  experimental  case  (unfortunately).  |<ymJ  <3  1 
and  the  energy  gain  Gm  on  the  mode  m  becomes 


f'm  =  2  Re  [i, . j  =  j’  l~^-  2»  £,•  I 

n  ur. 


-cosv.—  -,sinvt 


This  is  exactly  the  gain  one  calculates  by  using  the 
filling  lactor  obtained  by  integrating  the  mode  profile 
overlap  with  the  gain  profile.  Specializing  to  the 
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TEM00  case  with  a  Gaussian  electron  beam  of  width  <r, 
we  find 


complete  with  the  familiar  filling  factor. 

The  vc  dependence  of  Gm  is  the  well  known  spectral 
dependence.  The  imaginary  part  of  gmn  is  not  new.  It 
describes  the  phase  shift  of  the  radiation  field  as 
described  by  Colson  [18].  The  inhomogeneous 
broadening  term  in  (31)  clearly  distorts  and  reduces  the 
magnitude  of  the  gain  spectrum  if  it  is  present  in  the 
integral  of  (30). 

The  effect  of  the  divergence  of  the  beam  on  the 
diagonal  terms  in  G  is,  to  first  order ,  and  for  a 
filamentary  electron  beam,  the  addition  of  a  time- 
varying  phase  which  shifts  the  resonance  curve  in  (32) 
by  a  constant  depending  on  the  mode 

V‘^V‘“^(2p  +  ,+  1)-  (35) 
7tW0 

Equation  (35)  means  that  the  gain  curves  of  the  modes 
are  shifted  with  respect  to  each  other.  This  effect  has 
been  calculated  for  the  fundamental  TEM00  mode  in 
the  energy  loss  approximation  [14],  and  has  recently 
been  observed  experimentally  at  Orsay  [21].  It  should 
be  noted  that  for  many  practical  situations  where  the 
cavity  is  optimized  for  gain  on  the  TEM00  mode,  this 
expression  is  valid  for  only  the  lowest-order  mode.  The 
higher  modes  become  distorted  in  form  as  well  as 
simply  shifted  in  resonance  parameter  by  (35). 


2.2.  The  Low-Field  Oscillator 

We  now  discuss  some  consequence  of  the  linearity  of 
the  low-field  problem  on  the  optimization  of  an  optical 
cavity  for  an  FEL  oscillator  experiment.  In  such  an 
experiment  the  light  pulses  reflect  n  times  on  the  cavity 
mirrors  («^2)  between  interactions  with  electrons  in 
the  undulator.  The  matrix  governing  the  mode  evolu¬ 
tion  from  one  amplification  to  the  next  is 

(/  +  G)C,  (36) 

where  G  in  the  gain  matrix  defined  previously,  and  C  is 
the  cavity  matrix  describing  the  n  reflections  on  the 
mirrors.  In  a  set  of  cavity  eigenmodes.  C  is  diagonal 

Ctt'-AV v\  (37) 

with  cjj  =  n1  =^exp(ix;),  where  1  - gj  are  the  total 
losses  on  the  »:  reflections,  including  transmission, 
absorption,  scattering,  and  diffraction.  If  diffrac¬ 
tion  is  negligible,  the  eigenvectors  v1  become  the 


Gaussian  TEM,,,  modes,  and  the  phase  shift  per  round 
trip  x,  becomes,  for  n  =  2  reflections  per  amplifica¬ 
tion  and  identical  radius  of  curvature  mirrors. 
x(1,  =  4(2p  +  /+  l)tan  2r„).  where  /.,  is  the  optical 
cavity  length. 

The  matrix  (36)  is  the  fundamental  matrix  of  the 
problem.  Us  diagonalization  allows  the  calculation  of 
the  mode  evolution  up  to  the  onset  of  saturation: 

(/  +  G)C  =  PAP~  1 ,  (38) 

[(/  +  G)C]m  =  P/PP'1,  (39) 

where  the  columns  of  P  are  composed  of  the  eigenvec¬ 
tors  of  (I  +  G)C,  and  ,1  is  diagonal.  The  fastest  mode- 
growth  will  be  obtained  with  the  eigenmode  having  the 
highest  eigenvalue  modulus.  Optimization  of  the  FEL 
oscillator  will  then  consist  of  maximizing  the  desired 
eigenvalue  of  U  +  G)C. 

If  the  gain  is  low  (as  it  is,  unfortunately,  for  our  system 
on  ACO),  one  can  diagonalize  the  evolution  matrix 
(36) 

(I +  G)C:‘ =  (40) 

using  the  cavity  eigenmodes  x‘  as  the  basis  for  a 
perturbation  expansion  of  the  new  modes  To  first 
order  in  the  non-degenerate  case,  the  result  is 


~  -  je'Ij(  1  +  Qjj)  • 


-j  _  yj  y  _ — m _ _ _ (.  Y" 


(41) 

(42) 


Under  these  conditions,  the  FEL  design  is  optimized 
by  maximizing  the  diagonal  term  gJt  corresponding  to 
the  desired  mode.  From  (33)  and  (12)  it  is  clear  that  the 
beam  size  iv0  must  be  reduced  down  to  the  order  of  the 
electron-beam  size  in  order  to  optimize  the  coupling, 
but  if  the  mode  becomes  too  divergent,  the  time 
dependent  terms  in  Em  and  £„  of  (28)  begin  to  reduce 
the  gain.  The  optimal  situation  lies  between  these  two 
extremes,  and  has  been  calculated  in  detail  (using  the 
energy  loss  approximation)  by  Colson  and  Elleaume 
[14]. 

The  optimization  procedure  must  also  be  limited  by 
the  stability  condition  [17]  on  the  cavity.  For  the 
Orsay  experiment,  the  radius  of  curvature  chosen  to 
optimize  the  small-signal  gain  was  P  =  3  m.  which  is 
acceptably  close  to  the  stability  limit  of  2.75  m.  There 
are  cavity  designs  in  which  C  is  degenerate  for  which 
the  optimization  procedure  is  not  necessary.  For  these 
designs,  Xj  is  a  constant  independent  of  the  index.  In 
these  cavities,  any  combination  of  modes  reproduces 
itself  after  n  reflections.  If  n  =  2  as  in  the  Stanford  and 
the  Orsay  experiments,  the  concentric  and  the  plane- 
parallel  cavities  are  degenerate,  and  the  eonfocal  cav  ity 
is  degenerate  on  the  p  modes  (quasi-degenerate).  These 
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Fig.  I.  Simplified  schematic  diagram  of  the  gain-measurement  ap¬ 
paratus  [23]  showing  the  argon-laser  focussing  system,  the  collimat¬ 
ing  iris,  and  the  double  demodulation  detection  system 


cavity  designs,  however,  are  useless  since  they  stand 
critically  on  the  stability  boundary.  The  tolerance  on 
the  mirror  radius  of  curvature  is  on  the  order  of  |t/00[ 
(obtained  from  the  perturbation  expansion)  which  is 
difficult  to  meet  if  the  gain  is  low.  For  two-mirror 
devices  where  it  >2.  such  as  the  Novosibirsk  experi¬ 
ment  where  n  =  8.  in  general  for  mirrors  of  equal  radius 
of  curvature  there  exist  n,  2+1  cavity  designs  with  a 
degenerate  C  matrix: 


_  nm 
2  tan  — 

it 


(43) 


and  ii/2  with  a  quasi-degenerate  C  matrix  in  which  the 
odd  I  modes  change  sign  on  every  amplification.  Only 
two  of  the  degenerate  and  one  of  the  quasi-degenerate 
cavities  correspond  to  the  unstable  cavities:  the  others 
are  potentially  useable  in  an  experiment.  The  value  of  a 
degenerate  C  matrix  is  that  the  eigenvectors  of  the 
amplifier  plus  cavity  matrix  (36)  are  equal  to  those  of 
the  gain  matrix  alone,  multiplied  by  a  constant.  This 
degeneracy  allows  the  cavity  to  oscillate  on  the  most 
favorable  combination  of  modes  which  best  fits  the 
electron  beam  shape.  In  this  manner,  the  gain  can  be 
increased  by  factors  of  two  or  three  over  the  gain  of  an 
optimized  TEM„n  mode,  particularly  if  the  electron 
beam  size  is  smaller  than  the  TEM„„  mode.  The 
tolerance  on  the  mirror  radius  for  the  deeeneraev  of  C 


Fig.  2.  The  measured  gain  as  a  function  of  the  iris  diameter  [16] 
normalized  to  the  measured  beam  waist  at  the  iris.  The  solid  points 
were  taken  closing  the  iris  and  the  open  points  while  opening  n  The 
error  bars  are  the  one  sigma  statistical  errors  All  points  have  the 
same  horizontal  error  bar  which  is  shown  for  the  point  at  2.7.  The 
solid  curve  is  calculated  using  the  measured  values  for  the  electron 
and  laser  beam  sizes.  The  effect  of  each  higher-order  mode  is  show  n 
by  the  dashed  curves 


will  still  be  tight,  and  the  experimental  utility  of  these 
cavities  remains  to  be  invesliaated. 


3.  Application  to  Gain-vs-Apcrturc  Experiment 

3.1.  Description  of  the  Experiment 


The  gain  of  the  Orsay  FEL  has  recently  been  measured 
with  the  optical  klystron  in  place  [22]  in  an  amplifier 
experiment  using  an  external  argon  ion  laser  to  pro¬ 
vide  the  coherent  mode.  A  detailed  description  of  the 
apparatus  can  be  found  in  [23],  and  a  schematic 
description  is  given  in  Fig.  1.  The  laser  beam  is  anal¬ 
yzed  at  a  distance  </  from  the  optical  klystron  after 
passing  through  an  adjustable  collimating  iris  (Fig.  1), 
which  is  centered  on  the  laser  mode  emerging  from  the 
interaction  region.  The  gain  is  measured  as  the  ratio  of 
the  power  detected  in  phase  with  both  the  electron 
repetition  frequency  and  the  chopper  frequency  (the 
amplified  power)  divided  by  the  power  in  phase  with 
the  chopper  alone  (the  incident  laser  power). 
Calibration  is  performed  as  in  [23], 

The  again  is  recorded  as  a  function  of  the  iris  aperture, 
and  large  variations  are  observed  [  16].  One  set  of  data 
points  is  reproduced  in  Fig.  2.  where  the  gam  is 
normalized  to  its  value  for  the  iris  completely  open, 
and  the  iris  diameter  is  normalized  to  the  measured 
beam  waist  at  the  iris.  The  data  is  taken  at  maximum 
gain,  which  means  r  id  for  the  optical  klystron,  and 
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Fig.  3.  Calculated  curves  for  the  gain  as  a  function  of  iris  diameter 
under  the  conditions  of  the  Orsay  experiment  [16],  The  electron 
beam  dimension  I  -  a  (/ it/XL  is  varied  to  show  the  effects  of  the 
beam  size  on  the  excitation  of  the  higher  order  modes.  The  value 
I  =0.76  corresponds  to  a  =  0.35  mm  which  is  very  close  to  the  value 
at  which  the  experimental  points  were  recorded 
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Fig.  4.  Calculated  gain  as  a  function  of  iris  diameter  for  several  iris 
positions  d.  under  the  conditions  of  the  Orsay  experiment  [16],  The 
ratio  of  the  iris  to  optical  klystron  distance  J  divided  by  the  optical 
klystron  length  L  is  varied  through  the  range  05  to  x.  The 
experimental  points  of  Fig.  2  were  taken  for  dL  =  9 


the  laser  beam  was  carefully  aligned  to  within  about 
0.05  mm  of  the  axis  of  the  electron  beam.  The  change  in 
the  measured  gain  as  the  iris  is  closed  means  that  the 
laser  is  not  uniformly  amplified  in  its  transverse  profile. 
In  fact,  this  experiment  provides  a  very  sensitive 
technique  for  measuring  the  power  emitted  into  the 
higher-order  modes  even  in  the  small  gain  limit  and  for 
a  monomode  input  beam.  Clearly  a  calculation  of  the 
gmm  is  necessary  in  order  to  explain  these  results.  In  the 
next  section,  we  apply  the  theory  we  have  developed  to 
the  case  at  hand,  and  in  Sect.  3.3,  precise  comparison  is 
made  between  the  experimental  and  the  theoretical 
results. 


3.2.  Multimode  Emission  in  a  Single- Mode 
Amplifier  Experiment 

In  this  subsection,  we  assume  the  incident  wave  is  a 
single  mode  TEM00  beam  with  a  weak  field  (|u0|<  1). 
and  perfectly  aligned  onto  the  electron  beam.  As 
discussed  previously,  we  take  the  cylindrical  eigen- 
modes  based  on  the  form  of  the  input  beam.  Using 
the  notation  of  Sect.  1,  the  input  laser  field  reads 

£'(r)  =  c0£0(r)eiVo,,\  (44) 

where  the  subscript  0  refers  to  the  TEMno  mode  of  (12) 


and  (13).  From  (29),  the  output  field  £’(r)  becomes 

£5  =  £'+f0  £  gj0EJ(r)e‘'-^.  (45) 

j-o 

Assuming  low  gain,  the  output  power  passing  through 
the  iris  aperture  is 


—  =  J  </S|£!|2  *  cl  i  dSEl  +  2 cl  Re 
•[|o9Joi^£0£Jc'---j, 


where  JdS  covers  the  iris  aperture.  The  gain  is 
therefore 


2  Re  I  (/,J</S£0E,e"--v 
G  = - ^2 - = - — - 

J  ‘ISEl 


For  purely  cylindrical  1  =  0  modes.  G  can  be  written 


G  =  2  Re  j  (jno  +  V  g  e 


:  )  Ux)cxdx 

t  2/>i  jn  1  —  Q 

e  -o - 

P"  V 


where  is  the  Rayleigh  range  of  the  laser  mode,  r  is 
the  distance  between  the  ins  and  the  laser  beam  waist, 
and  ,V  =rl  2iv:(r)  where  r0  is  the  iris  diameter.  There 
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are  two  interesting  limiting  cases: 

G  =  2  Re { £/00 }  X-+  tc  (iris  open).  (49) 

G  =  2Rej(c/00)  +  £  gp0e'PUn  :o| 

X  -*0  (iris  closed) .  (50) 

It  is  obvious  from  (48-50)  that  the  gain  changes  with 
iris  diameter  in  a  way  which  depends  on  the  magni¬ 
tudes  of  the  off-diagonal  terms  in  the  gain  matrix. 

The  generalization  is  straightforward  to  the  case  of  the 
multimode  input  beam,  and  to  imperfect  alignment  of 
the  laser  and  electron  beams,  although  the  calculation 
becomes  more  difficult.  This  calculation  also  applies  to 
high-power  input  laser  beams  (a0>l)  provided  one 
keeps  in  mind  that  the  gp0  are  functions  of  a0. 

3.3.  Application  to  the  Orsay  Experiment 

The  experimental  points  shown  on  Fig.  2  were  taken 
under  the  following  approximate  conditions 

laser 

beam:  -  measured  beam  waist  n0  =0.67 mm 

-  wavelength  =  5 145  A 

-  measured  beam  waist  at  iris  u(c)  =  2.7  mm 

-  distance  from  optical  klystron  to  iris 
d-  11.6  m, 

electron 

beam :  -  Gaussian  and  cylindrical  with 

<7=:  0.32  mm, 

optical 

klystron:  -  Nd  =  80  [19,  22.  24] 

-  resonance  parameter  corresponding  to 
maximum  gain  with  iris  open. 

The  solid  curve  of  Fig.  2  has  been  calculated  using  (27, 
28.  and  50)  for  the  planar  configuration  (24  and  25). 
taking  into  account  the  10  lowest  order  1  =  0  modes. 
The  dashed  curves  of  Fig.  2  show  the  contribution  of 
each  individual  mode.  These  curves  are  the  same 
whether  an  undulator  or  optical  klystron  is  used.  Very 
similar  curves  (not  shown)  were  calculated  for  other 
resonance  parameters  indicating  that  as  expected,  the 
diffraction  effects  do  not  change  much  as  a  function  of 
detuning  parameter  for  modes  with  low  divergence. 
Figure  3  shows  the  calculated  effect  for  several  dimen¬ 
sionless  electron  beam  transverse  dimensions 
2T  =  a n!kL  =0.4.  0.76.  1.5.  3.  where  L  is  the  length  of 
the  magnetic  interaction  region.  1=0.76  corresponds 
to  the  value  o  =0.35  mm,  close  to  that  used  in  Fig.  2. 
The  flattening  of  the  curves  as  2,'  is  increased  to  E  =  3  is 
due  to  the  vanishing  of  <yn, ‘/no  f°r  7*0  as  a 
increases. 


Figure  4  shows  the  calculated  effect  for  various  iris 
distances  d  from  the  optical  klystron  center,  normal¬ 
ized  to  the  optical  klystron  length :  J  L  =0.5.  1.  1.8,  2.5. 
and  jc.  The  experimental  points  of  Fig.  2  were  obtain¬ 
ed  for  d/L  =  9.  The  inversion  of  the  effect  is  due 
primarily  to  the  term  exp[i2ptan~  '(c/r0)]  with  p=  1  in 
(50)  which  switches  from  -t- 1  to  —  1  as  r  goes  from  zero 
to  infinity  (the  mode  p=  1  gives  the  predominant 
effect).  At  short  distances  the  TEM,0  mode  interferes 
constructively  on  axis  1,  and  at  long  distances,  it 
changes  sign. 
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ABSTRACT 


Me  report  the  first  measurement  of  the  off-diaaonal  terms  in 
the  transverse  aain  matrix  in  a  free  electron  laser.  The  hiaher 
order  transverse  modes  stimulated  in  the  FEL  interaction  are  shown 
to  diminish  in  amplitude  as  the  electron  beam  size  is  increased, 
and  to  be  insensitive  to  the  resonance  parameter.  Remarkably  close 
agreement  is  demonstrated  with  the  theory  Cl  3.  The  effects  of  a 
multiple  mode  input  beam  are  also  measured,  and  the  theoretical 
p.,pression  is  derived. 


INTRODUCTION 


Besinnina  in  1980  a  series  of  measurements  have  been 
performed  of  the  stimulated  amplification  of  the  Orsay  storaae 
rina  free  electron  laser.  These  measurements  were  part  of  the 
successful  effort  to  achieve  oscillation  us i ns  the  ACO  storaae 
rina  system  C2r33.  The  original  aain  experiments  14.53  on  the 
superconducting  undulatorf  and  a  later  series  of  measurements 
CG.73  on  the  permanent  maanet  undulator  NOEL  and  its  optical 
Klystron  conf isurat i on  demonstrated  that  the  amplitude  and  the 
spectrum  of  the  aain  could  be  correctly  calculated  by  the 
classical  independent  particle  theory  C83.  These  experiments  Merc 
sensitive  to  the  total  energy  emitted  by  the  stimulated  emission 
rrocessr  and  did  not  provide  any  information  on  the  spatial 
structure  of  the  amp  1 i f i cat i on  process.  In  this  paper  we  report  on 
an  extension  of  these  measurements  which  probes  the  three 
dimensional  nature  of  the  interaction. 

By  means  of  an  adjustable  aperture  placed  in  the  amplified 
beam.'  we  have  been  able  to  measure  the  excitation  of  hiaher  order 
|  transverse  modes  induced  by  the  nonuniform  transverse  profile  of 

the  electron  beam.  This  information  is  mapped  into  the  spatial 
distribution  of  the  electric  field  at  the  analyzing  aperture.  The 
|  ratio  of  the  transmitted  stimulated  power  to  the  transmitted  power 

I 

t 


of  the  external  laser  is  measured  as  a  function  of  the  aperture 


diameter.  The  functional  dependence  of  the  ratio  on  the  aperture 
diameter  1:  a  diagnostic  on  the  transferee  mode  mixing:  which  has 
occurred  durina  the  amplification  process. 

If  the  spatial  distribution  of  the  stimulated  radiation  is 
identical  to  that  of  the  incident  laser?  their  diffraction 
characteristics  are  identical?  and  the  ratio  of  the  powers 
transmitted  throuah  the  aperture  will  be  independent  of  its 
diameter.  A  dependence  will  be  seen  only  if  the  stimulated 
radiation  contains  a  different  combination  of  transverse  modes 
than  the  incident  radiation.  The  aperture  dependence  of  this  ratio 
is  therefore  a  ssood  rfiaanottic  of  transverse  mode  mixine  in  the 
laser  amplifier. 

In  a  free  electron  laser?  the  input  laser  mode  is  perfectly 
renlicated  only  under  the  conditions  that  the  electric  field 
strenath  remains  below  the  saturated  level?  that  the  electron  beam 
is  uniform  across  the  optical  mode?  and  that  the  divergence  of  the 
laser  is  nealiaible  in  the  interaction  rosiori.  If  the  electron 
beam  dimensions  are  comparable  to  those  of  the  optical  mode?  if 
saturation  occurs?  or  if  diffraction  is  significant?  higher  order 
modes  will  be  excited.  In  most  configurations  of  practical 
interest?  one  or  most  of  the  latter  conditions  will  hold. 

The  modo  content,  of  the  optical  field  propagating  through  the 
PEL  van  bo  conveniently  thought  of  as  a  vector  whose  elements  are 


the  amplitudes  and  phases  of  the  modes.  As  the  beam  pro  pa sates 
throuah  free  space-  the  phases  evolve  in  a  fashion  characteristic 
of  each  mode?  but  the  amplitudes  remain  constant  <in  the  absence 
of  lossesl.  The  amplifier  contributes  an  additional  vector  of  complex 
field  amplitudes  which  describes  the  stimulated  radiation/  and 
which  is  in  aen  era  1  not  a  simple  multiple  of  the  input  vector  C 1  1  . 

The  Field  vector  at  the  output  of  the  FEL  amplifier  is  the  sum  of 
the  input  vector  propasated  to  the  output  and  the  stimulated 
emission  vector.  The  maanitude  squared  of  the  output  vector  will 
be  I arser  than  that  of  the  input  by  a  factor  which  is  the  net 
amplification/  and  the  vector  will  be  rotated  in  the 
multi-dimensional  mode  space.  If  the  sain  is  larso  and  mode  m  i  x  i  na 
occurs/  the  output  field  vector  will  be  stronsl '/  rotated.  Under 
these  circumstances.-  the  laser  can  be  expected  to  run  multi-mode/ 
and  a  self-consistent  calculation  of  the  field  distribution  must 
be  performed  to  identify  the  steady  state  field  vector  and  the 
a a  in  Jt  experiences.  In  a  low  sain  oscillator/  however/  the  same 
level  of  mode  mixina  will  produce  no  effects  because  the  rotation 
of  the  vector  is  nealiaible. 

The  aai  n  of  the  G  r  s  a  */  FEL  is  very  loWr  on  the  order  of  luA-3 
C -  7 1 .  Even  thoush  the  conditions  for  mode  mix  ins  are  present  (the 
dimensions  of  laser  and  electron  beam  are  comparable).-  the 
rotation  of  the  field  vector  is  minuscule.  As  could  be  expected 
under  those  conditions  r  the  oscillator  produced  a  nearly  pure 
TEMoo  Gaussian  mode  C31.  (Even  the  Stanford  oscillator  which  has  a 


sain  of  about  10%  C93  produces  a  Gaussian  mode  CIO]).  The  newer 
FEL  systems  now  under  construction  at  "Etanf ordr  Orsa1/?  and 
elsewhere  are  projected  to  have  aains  per  pass  ransms  aboue 
unit:*.  The  mii'ina  phenomena  of  the  transverse  modes  will 
play  an  imrortant  role  in  the  operation  of  these  new  devices.  Our 
measurement  technique  permits  us  to  obtain  some  information  on  the 
three  dimensional  ampl if i cat i on  problem  even  in  devices  where  the 
sain  is  so  low  that  no  effects  can  be  observed  in  the  modes  of  the 
oscillator. 

Two  series  of  experiments  have  been  performed?  both  with  the 
undulator  NOEL  modified  into  the  optical  Klystron  configuration 
C7J.  The  first  series  of  experiments  (which  we  will  identify  as 
series  J>  served  to  identify  the  phenomenon  ?  and  stimulated  our 
theoretical  effort  Til  on  the  problem.  A  second  series  (series  II) 
was  planned  in  order  to  extend  the  comparison  between  theory  and 
experiment  b  •/  chanains  the  placement  of  the  analyzing  aperture.  We 
report  here  on  the  results  of  these  two  experiments. 
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EXPERIMENTAL  METHOD 


The  apparatus  we  used  was  essentially  identical  to  that  reported 
in  the  original  sain  measurements  C4-71.  Those  unfamiliar  with  the 
details  of  the  setup  should  refer  to  the  schematic  diaaram  in  T43 
alons  with  the  associated  description  of  the  experiment  and  the 
calibration  proceed u re.  A  chopped  external  arson  laser  is 
focussed  coaxially  onto  the  electron  beam  in  the  optical  klystron 
where  it  is  amplified  by  the  FEL  interaction.  The  output  beam  is 
passed  through  an  adJustible  aperture.-  filtered  in  a 
mo nochroma tor?  and  detected  in  a  resonant  detector.  The  amplified 
power  is  extracted  in  a  double  demodulation  schemer  a  12  Mhz 
lock-in  detects  the  sum  of  the  spontaneous  and  the  stimulated 
power.-  and  a  subsequent  lock-in  opera',  ins  at  the  frequency  of  the 
chopper  extracts  a  si sinal  proportional  to  the  stimulated  power  at 
the  detector. 

In  this  experiment:  the  undulator  sap  is  scanned  to  set  a 
maximum  of  the  optical  klystron  33in  spectrum  C7D  at  the  laser 
inuelensth.  Alianment  is  performed  by  ortimirina  the  amplitude  of 
the  stimulated  power  when  adjusting  the  transferee  angles  and 
positions  of  the  electron  beam.  A  half-wane  plate  in  the  laser 
input  beam  is  also  adjusted  to  ensure  th3t  the  polarization  of  the 


input  mode  is  parallel  to  that  of  the  undulator  emission. 


One  anslvsins  aperture  was  set  up  as  close  as  possible  to  the 
interaction  resjonl  Just  outside  the  e x  i  t  window  at  a  distance  o F 
288  cm  from  the  center  of  the  optical  Klystron.  This  unit  was  a 
set  of  fiyed  apertures  of  diameters  ?  I.O.-  1.5?  2.0r  and  2.5  mm/ 
mounted  on  a  rotating  exchanger  blocK.  The  selected  aperture  could 
bo  rotated  into  position  and  aligned  precisely  in  the  plane 
transverse  to  the  laser.  Before  each  measurement  with  this  set  of 
apertures.-  alignment  was  established  by  reeenterina  the  diffracted 
beam  onto  the  original  axis.  A  second  aperture  was  placed  far  from 
the  interaction  region  at  a  distance  of  14.7  m.  Since  the  laser 
beam  was  larger  heror  a  minimum  aperture  of  2  mm  was  sufficient/ 
and  an  adjustable  iris  was  used:  eliminating  the  need  for 
realignment  on  each  change  in  diameter.  The  powers  of  the 
transmitted  laser  arid  of  the  stimulated  omission  were  measured  in 
separate  channels  and  recorded  for  about  100  seconds.  The  laser 
was  then  blocKed  in  order  to  permit  a  measurement  of  the  noise 
bacKaround?  the  aperture  diameter  was  readjusted?  the  aperture  was 
realigned  (if  necessary)-  and  the  two  power  measurements  were 
repeated.  The  curves  of  figures  4-8  were  taken  by  this  method. 

In  order  to  verify  that  no  detection  artifacts  were 
masquerading  as  real  effects?  we  verified  first  the  linearity  of 
the  detector.-  and  then  the  position  sensitivity  of  its  high  and 
low  frequency  response.  B“  inserting  a  series  of  broadband  Filters 
up  to  a  neutral  density  of  we  simulated  the  reduction  in  power 


produced  by  closing  the  aperture  diameter  without  chansing  the 


r^Uo  of  laser  to  stimulated  power.  The  measured  ratio  remained 
within  a  fraction  of  one  standard  deviation  from  unity.  The 
detector  was  also  mooed  from  side  to  side  with  a  focussed  input  to 
checK  that  spatial  variations  in  its  sensitivity  were  not 
responsible  for  the  sisnal.  The  same  nesative  result  was  obtained! 
the  detection  system  is  linear  and  sensitive  only  to  the  power 
transmitted  throuah  the  analyzins  iris. 

Since  our  soa!  is  to  measure  the  dependence  of  the  stimulated 
Power  on  the  aperture:  it  i?  particularly  important  to  be  sure 
that  no  portion  of  the  laser  output  mode  is  scraped  off  at  any 
point  in  the  optical  transport  system.  'The  aperture  derendance 
disappeared  when  the  monochromator  slits  were  closed  down: 
scrapina  o *■' f  a  portion  of  the  ^ocal  area).  Durina  set-uP:  the 
laser  mode  was  centered  on  the  transport  mirrors r  and  the 
monochromator  slits  were  opened  to  2  mm  to  ensure  that  when  the 
analyzins  operture  was  open>  all  of  the  stimulated  power  entered 
the  detector.  Mhon  the  analyzins  aperture  was  placed  in  its 
farthest  position  From  the  undo  1 ator >  no  problem  was  observed  in 
power  transport.  Mhen  the  aperture  was  placed  in  its  near  position 
to  the  undulatorr  and  closed  to  its  minimum:  the  srowth  of  the 
output  beam  duo  to  diffraction  did  result  in  perhaps  a  10*4  power 
loss  on  one  of  tho  beam  transport  mirrors. 

The  interpretation  of  the  results  demands  an  accurate 
fnouledse  of  the  electron  and  optical  beam  sizes.  In  an  attempt  to 


simplify  the  comparison  with  theory?  considerable  effort  was 
expended  to  make  the  external  Arson  laser  run  in  a  single 
transverse  mode.  For  series  I?  this  effort  was  successful.  For  the 
later  series?  the  asod  laser  plasma  tube  performed  less  well?  and 
even  with  an  iris  internal  to  the  Arson  laser  oscillator  closed 
and  optimized?  the  beam  showed  significant  higher  order  mode 
rower . 

The  laser  beam  sizes  wore  repeatedly  measured  before  and 
during  series  II.  The  input  laser  mode  was  measured  with  a  Reticon 
diode  array  with  a  25  micron  diode  separation.  The  beam  shape  and 
the  full  width  at  half  maximum  was  recorded  at  five  positions 
about  the  focus?  both  in  the  vertical  and  the  horizontal 
dimensions.  Those  measurements  were  fit  to  a  Gaussian  TPMoo  mode 
to  extract  the  beam  waist  and  focal  position  information.  The 
fits?  shown  in  figure  1?  show  moderate  agreement  with  the  data? 
but  cannot  be  used  reliably  to  estimate  the  higher  mode 
intensities.  The  astigmatism  of  the  Arson  laser's  asymmetric  cavity 
is  apparent  •  as  expected.-  the  vertical  focus  is  tighter  than  the 
horizontal?  and  occurs  slightly  before  it.  Figure  2  shows  the 
transverse  profile  of  tho  intensity  at  the  center  of  the  undulator 
(nearly  Gaussian)?  and  4.3  m  before  the  focus?  where  the  higher 
mode  power  becomes  apparent. 

For  the  data  analysis?  we  also  require  a  Knowledge  of  the 
size  of  the  laser  beam  at  the  aperture.  This  information  is  obtained 


from  the  dependence  of  the  total  detected  laser  power  on  the 
diameter.  If  the  diameter  is  Known  and  the  aperture  is  well 
aligned:  the  laser  beam  size  can  be  calculated  ‘■'rom  the  reduction 
of  the  transmitted  laser  power:  under  the  assumption  of  a  Gaussian 
boam. 


The  electron  beam  size  is  less  well  Known.  It  is  measured  as 
shown  in  figure  1.  The  synchrotron  radiation  emitted  from  the  beam 
in  the  dipole  magnet  15S  cm  in  front  of  the  center  of  the  optical 
Klystron  is  imaged  onto  the  Reticon.  The  vertical  and  the 
horizontal  boam  shapes  are  typically  observed  to  be  Gaussian.-  with  a 
width  which  "aries  according  to  the  current  in  the  ring.  These 
dimensions  also  change  with  the  tunes  of  the  ring.-  so  separate 
electron  beam  measurements  are  required  on  each  injection  during 
the  oy per 3 merit  .  Ho  diagnostics  are  available  on  the  beam  size  in 
the  center  of  the  optical  Klystron.  This  is  an  unfortunate 
limitation  since  the  beta  functions  of  ACO  have  not  been  measured 
with  the  soxtupoles  in  operation.  In  the  absence  of  the  trajectory 
errors  induced  by  the  sextupolesr  the  theoretical  model  of  the  ACO 
ortics  predicts  that  the  beam  sizes  change  by  less  than  25V.  and 
6\r  respectively  in  the  horizontal  and  the  vertical  planes.-  when 
moving  from  the  electron  beam  measurement  point  to  the  center  of 
the  OK.  These  numbers  can  also  be  taKen  as  a  rough  estimate  of  the 
systematic  error  in  the  electron  boam  size  determination. 


Due  to  a  technical  problem  with  the  extraction  mirrorr  the 


second  series  oF  experiments  was  performed  with  the  probe  laser 
beam  exiting  the  vacuum  system  throush  the  rear  oscillator  mirror. 
Tho  transmission  of  this  mirror;  a  hish  reflector  at  6228  A.-  is 
about  80/i  at  the  Arson  51^5  A  wauelensth.  It  acts  as  a  defocussins 
lens  with  a  focal  lensth  of  about  -5.2  m.  These  experiments  also 
had  to  bo  performed  at  low  enerax  '!166  Met*)  in  order  to  avoid 
damasins  C 1 1 1  the  oscillator  mirrors  with  the  L,,.,  senerated  in  the 
optical  Klystron.  All  measurements  were  done  with  the  5145  A  laser 
line"  in  Series  I.-  Nd  =  80.-  K  -  1  .89.-  and  Wo  =  670  microns*  in 
Series  II-  Md  =  70;  K  =  .9;  and  Wo  =  610  microns. 


RESULTS  OF  THE  EXPERIMENTS 


The  most  precise  data  scan  taP.en  in  series  I  is  illustrated 
and  compared  with  theory  in  fisure  4.  Here:  the  input  laser  beam 
is  a  pure  TEMoo  mode:  and  the  amplified  beam  propaaates  freely  for 
14.3  m  before  being  analyzed  with  the  iris.  The  error  bars  can  be 
estimated  from  the  accuracy  with  which  the  two  consecutive  scans 
overlap.  The  solid  curve  is  the  theoretical  calculation  C13  based 
on  tho  measured  distances  and  beam  sizes?  no  fittina  has  been 
done.  Tho  agreement  is  excellent. 


This  is  the  first  experiment  to  demonstrate  the  existence  of 
the  off-diasonal  terms  in  the  gain  matrix  (transverse  mode 
mixing)?  and  establishes  the  validity  of  our  theoretical 
under; tandins  of  mode  mixing  in  FELs  in  the  weak  field  limit.  From 
the  theoretical  curve?  the  amplitudes  of  the  first  few  terms  in 
the  matrix  are  determined  to  be  \*A-l  =  .46?  and  P’/V.|  =  .17? 
the  hiaher  terms  are  less  than  10%.  It  is  clear  that  sianificant 
mode  mixina  is  ocurrina  with  an  electron  beam  to  laser  beam  size 
ratio  of  ^  ~  2<T/Mo  =  0.96.  If  a  laser  of  moderately  hish  aain  is 
crerated  with  such  a  beam  ratio?  the  oscillator  mode  will  differ 
significantly  from  a  Gaussian. 


Figure  5  show  the  results  of  two  series  I  measurements  taken 
as  tho  stored  current  decayed  by  a  factor  of  1.9?  the  resonance 
parameter  was  changed  by  a  very  small  amount  ($x=.l)  between  the 


too  scans.  Sianificant  differences  are  apparent  between  the  shapes 
of  these  curves.  The  error  liars  can  be  estimated  from  the 
deviations  in  the  high  current  measurements?  for  this  scan.-  two 
complete  cycles  of  the  aperture  were  made. 

This  data  demonstrates  the  effect  of  a  reduction  of  the 
electron  beam  size.  In  ACQ  the  transverse  dimensions  are  stronaly 
dependent  on  the  current  stored  in  the  ring  because  of  the 
multiple  TouscheK  effect  and  ion  trappins.  The  electron  beam 
dimensions  have  fallen  from  approximately  <T  =  620  microns  to  <T  =  4 
microns  between  these  two  data  scans,  which  makes  the  beam  ratios 
1.9  and  1  .  * ,  respectively.  The  curve  with  the  larsost  current 
(79  ma)  shows  less  gain  depression  at  small  aperture  then  the  low 
current  curve.-  taken  at  42  ma.  The  conclusion  is  that  smaller 
electron  beam  sizes  result  in  more  excitation  of  the  higher  order 
transverse  modes.  The  shape  of  each  curve  depends  on  the  shape  of 
the  mode  integrals,  and  shows  some  variation  as  the  beam  size 
drops . 

Figure  6  shows  a  later  set  of  data  taken  after  the  lifetime 
of  the  stored  beam  had  stabilized  to  a  larger  value.  These  three 
sets  of  points  were  taken  consecutively  at  approx i ma t e 1 v  the  same 
electron  beam  size  (<T=  420?  260.-  and  350  microns),  but  widely 
varying  resonance  parameters  (V  =  3.G.  3.0  .and  -2.0).  which 
correspond  to  the  f  nurt  h-*-'rom-t  ho-center  gain  and  absorption  peaks 
on  either  side  of  the  optical  klystron  spectrum.  Ho  '-ystematic 


deviation  appears  between  those  cur"es. 


Th  i  s  conclusion  is  further  supported  b  ”  a  second  type  of 
measurement,  for  a  31  vert  aperture  diameter?  a  scan  was  taken  us. 
the  resonance  parameter.  The  peak  positions  were  identified?  and  a 
new  scan  was  taken  for  another  diameter.  No  spectral  shifts  or 
chanaes  in  the  shape  of  the  spectrum  were  observed  over  a  ranae  of 
iris  openings  between  full  open  and  .9  on  the  horizontal  scale  of 
figure  6. 

The  theory  C 1 3  of  the  mode  mixing  effect  describes  both  of 
these  phenomena.  The  amplitudes  of  the  off-diasonal  terms  in  the 
gain  matrix  (see  equation  <21  of  the  next  section)  depend  on  the 
transverse  integral  of  the  electron  beam  with  the  relevant  modes. 
All  off-diagonal  terms  drop  to  zero  in  the  large  electron  beam 
limit.  This  is  the  behavior  verified  in  figure  5.  It  is  also 
apparent  from  <21  that  in  the  low  diffraction  limit?  the  spectral 
term  has  separated  out  as  an  independent  factor.  This  implies  that 
the  relative  amplitudes  and  phases  of  the  gain  matrix  are 
independent  of  the  resonance  parameter.  Figure  6  supports  this 
resul  t . 

In  the  large  diffraction  case?  the  resonance  curves  are 
shifted  with  respect  to  each  other  Cl  I  due  to  the  divergence  of 
the  modes.  The  amplitude  of  the  cross  terms  in  the  sain  therefore 
depends  on  the  divergence  of  *ho  source  and  the  scattered  modes. 


f, 


p  o  r  our  case?  t h  e  divergence  of  the  first  feu  modes  is  v  e  r  v  1  o  i j , 
so  no  dependence  is  observed  on  the  resonance  parameter. 

The  data  of  series  II  were  taken  with  the  primary  aoal  of 
verifying  the  dependence  on  the  distance  of  the  aperture  from  the 
interaction  resion.  The  evidence  of  multiple  mode  input ?  however? 
severely  complicates  the  analysis  of  the  results.  Figure  7  shows  a 
composite  of  four  successive  measurement  scans  taken  at  a  distance 
of  2.8B  m?  only  13  cm  from  the  output  mirror.  The  larse 
fluctuations  in  the  amplitude  are  caused  by  the  need  to  realign 
the  aperture  each  time  a  new  opening  is  selected.  It  is  clear  that 
a  large  increase  in  the  gain  is  observed  at  small  aperture 
openings?  although  the  exact  shape  of  the  curve  is  not  well 
determined. 

The  solid  cur*<e  is  the  theoretical  calculation  based  on  the 
assumption  of  a  single  mode  input  with  the  measured  beam 
parameters.  The  electron  and  laser  beams  were  further  approximated 
to  be  round  although  both  were  slightly  elliptic?  as  apparent  from 
figure  1  and  the  data  in  Figure  7.  We  believe  the  observed  higher 
order  mode  content  of  the  input  beam  is  primarily  responsible  for 
the  deviation  of  the  theoretical  curve  from  the  experimental  one. 
The  large  deviation  of  the  gain  from  unity  at  small  aperture  can 
bo  explained  if  power  is  present  in  the  higher  modes  of 
the  input  beam.  Additional  interferences  appear  with  the  input 
moies?  and  extra  sets  of  stimulated  modes  are  produced-  one  for 


each  input  mode.  The  theoretical  treatment  (see  the  next  section) 
can  become  suite  complex.  In  the  absence  of  detailed 
information  on  the  input  mode  content  >  a  comparison  between  the 
thaor:*  and  the  experiment  can  bo  successful  only  if  one  happens  to 
choose  the  correct  input  mode  intensities  and  phases?  and  would  be 
of  limited  t* a  1  u e  . 

A  set  of  data  at  the  14.3  m  distance  was  also  taken  in  series 
Ilr  and  is  p resented  in  fissure  .  The  error  bars  for  this  set  are 
comparable  to  those  of  the  other  14.3  m  data  of  figure;  4-6.  Our 
results  at  this  distance  show  a  roughly  constant  or  increasing 
ratio  in  a  situation  whore  the  single  mode  theory  predicts  a  4  07. 
reduction  at  small  aperture  diameter.  The  two  scans  shown  were  taR 
at  different  electron  beam  sizes.-  and  again  show  the  reduction  of 
the  hi3her  order  terms  at  larger  beam  size  (see  figure  5). 

The  data  of  Figure  8  were  taken  under  the  same  conditions  as 
those  of  figure  71  only  the  electron  beam  sizes  are  slightly 
different.  Houeuer  there  is  a  large  and  significant  difference 
between  the  two  sets  of  measurements.  He  are  observing  the  effect 
that  we  set  out  to  see-  namely  that  the  relative  phases  of  the 
higher  order  modes  change  as  the:'  propagate  through  free  space.  If 
cur  interpretation  of  the  origin  of  the  dependence  of  the  gain  on 
the  diameter  is  correct-  a  change  in  the  distance  at  which 
the  aperture  is  placed  should  shift,  the  phase  of  the  coefficients 
equation  <G>.  and  change  the  shape  of  the  gain  vs.  iris  curve  in  a 
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MODE  MIMING  THEORY  APPLIED  TO  THE  GAIM  US  APERTURE  EXPERIMENT 


E '  I  e  a  u m  e  and  Deacon  Cl]  ha" e  identified  the  equations  of 
motion  for  the  excitation  of  the  hishor  order  t ran s verse  modes  in 


an  FELr  and  solved  them  in  the  low  field  case.  Their  principal  low 


field  results  are  that  the  complex  mode  amplitudes  evolve 
a c c o r d i n 3  to 


The  linear  relationship  between  input  and  output  mode  amplitudes 
holds  even  1  n  the  hish  sain  resion  Cl].-  althoush  the 

transformation  becomes  difficult  to  calculate.  In  what  follower  we 
restrict  our  examples  to  the  small  sain  limit.  In  this  limit  r  for 
low  diffraction :  and  for  small  beam  slippase 

parameter  (as  in  the  Orsax  experiment) r  the  elements  of  the  sain 
matrix  *3  are 
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Let  us  apply  these  equations,  folio wins  C 1 1 ,  to  the  case 
in  which  the  incident  beam  contains  an  arbitrary  number  of  modes. 
The  input  field  is 


£r  =  cj  £je 


A 


where  we  use  the  summation  convention:  repeated  indices  are 
summed.  The  functions  E<x.-y>»  and  ^  y.  ,  y )  are  defined  in  C  1 1 
as  the  normalized  amplitude  and  the  phase  of  the  mode.  The 
output  field  is,  lfrom  <  1  >  r 


C°  e 1  p  ' 

C  :  t  t  ESe  jj 


The  aain  measured  in  this  experiment,  which  is  the  ratio  of  the 
stimulated  power  to  the  incident  power  transmitted  t.hroush  the 
anerture,  is,  for  small  aain  s  y stems. 
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where  the  'ntearal  is  over  the  surface  of  the  aperture  opening. 


For  simplicity  of  notation,  let  us  define  a  normalized  mode 
integral 


which  defines  the  elements  of  a  Hermitian  matrix  T.  Since  the 
intearal  is  over  the  surface  out  to  the  aperture.-  it  is  clear  that 
these  quantities  depend  on  the  diameter.  Due  to  the  different  mode 
structures,  each  turn  has  its  own  aperture  dependence,  and  its  own 
unique  phase.  Substitutins  into  ( 6) the  sain  becomes 


This  relation  can  bo  expressed  more  compactly  in  matrix  notation. 
Adoptiris  C  as  the  vector  of  complex  mode  amplitudes,  we  find  the 
compact  result 

l  Ref  (cI(i+£)htg-ic)>  ] 


Physically  speaKina:  this  result  is  not.  hard  to  understand.  The 


numerator  is  the  stimulated  power  term.  Each  input  mode  Cj  produces  a 
set  of  stimulated  modes  uia  the  sain  matrix  ^  vector  of 

stimulated  modes  is  GIC>.  Each  stimulated  mode  "beats"  with 
another  mode  C ^  when  it  passes  throush  the  aperture?  this  other  mode 
is  either  an  incident  mode  <CII  or  an  amplified  mode  <C!G  . 

Each  combination  of  modes  produces  a  nonzero  contribution  to  the 
transmitted  power  when  the  aperture  diameter  is  finite.  This 
contribution  is  proportional  to  a  particular  mode  intesralr  whose 
amplitude  is  given  by  Tm  .  The  measured  sain  is  the  ratio  of  the 
stimulated  power  to  the  incident  power;  and  the  denominator  takes 
care  of  the  effects  of  the  iris  on  the  incident  beam. 

To  make  the  situation  more  clear/  let  us  look  at  the  explicit 
sums  in  the  case  of  two  input  modes r  "o"  and  "p".  He  assume  that 
the  input  beam  is  predominantly  "  o  "  »  and  linearize  in  Cp  .  Equation 
<S>  reduces  to 


(9) 


whore  wo  have  rationalized  the  denominator;  and  collected  the 
off— diasonal  terms  in  the  sums. 


Wo  note  here  that  for  sinsle  mode  input  <  Cp  =  0 > >  this  result 
reduces  to  that,  of  equation  (46)  in  C  1 1  .  The  additional  terms  are 
the  various  beating  terms  of  the  stimulated  emission  with  the 
input  nodes:  integrated  across  the  iris.  The  last  term-  for 
example:  is  the  product  of  the  mth  stimulated  mode  (produced  via 
<J«p  from  the  Pth  input  mode)  and  the  Oth  input  mode:  integrated 
across  the  analyzing  aperture-  and  normalized  to  the  input 
intensity.  It  is  the  last  throe  sums  which  are  apparently 
responsible  for  the  difference  between  the  experimental  points  in 
figures  7  and  8  and  the  single  input  mode  theory  of  C 1 3 . 


If  the  electron  beam  is  very  large  compared  to  the  laser 
beam:  the  off-diasonal  terms  <2>  of  G  go  to  zero  by  the 
orthogonality  of  the  modes.  The  diagonal  terms  all  have  the  same 
magnitude  due  to  the  mode  normalization.  Equation  (9)  reduces.- 
under  these  circumstances to 


G  = 


(Flat  electron  beam) 


(10) 


which  means  that  the  gain  for  a  combination  of  two  modes  is  the 

same  as  that  of  one  of  them  alone.  If  we  redefine  the  sum  of  two 

modes  as  cur  now  starting  mode:  it  is  clear  that  an  arbitrary 

third  mode  amplitude  can  be  added  without  changing  the  result.  It 
follows  by  induction  that  for  an  arbitrary  combination  of  modes.- 

the  direction  of  the  mode  vector  is  not  changed  by  the  FEL 


amplification  process  for  a  flat  beam 


A  measurement  of  the  sain  under  these  circumstances  will  show 
no  variation  as  a  function  of  the  transverse  position  in  the  beam 
'as  shown*  (10)  is  independent  of  iris  diameter).  Any  variation  of 
the  sain  in  the  transverse  plane  is  proof  of  mode  mijtins.  Althoush 
the  dat,3  of  fisures  7  and  8  were  taken  with  an  unknown  multiple 
mode  input*  the  deviation  from  unity  of  the  results  are  a  clear 
indication  of  the  presence  of  mode  mixins.  The  finite  beam  size 
allows  off-diasonal  terms  in  (2),  and  thereby  produces  the  mode 
mixing  observed  in  this  experiment. 


CONCLUSION 


This  experiment  is  the  first  capable  of  resold  ins  an 
excitation  of  the  hishor  order  transverse  modes  in  a  free  electron 
laser:  and  has  established  the  existence  of  off-diasonal  terms  in 
the  transverse  sain  matrix.  The  technique  we  have  developed  is  a 
useful  diagnostic  of  the  transverse  mode  mixing  effect  which  is 
invisible  in  today's  FEL  oscillators  due  to  the  low  sain :  but. 
which  will  i»lax  a  dominant  role  in  the  operation  of  the  high  gain 
devices  now  being  constructed  or  planned. 

High  gain  free  electron  lasers  will  shout  strong  effects  of 
mnde  mixing  in  anx  system  in  which  the  electron  and  laser  beam 
sizes  are  comparable.  We  have  demonstrated  that  mode  mixing 
effects  are  typically  large  in  free  electron  lasers  but  that  they 
have  been  mashed  in  the  existing  devices  by  the  low  value  of  the 
sain.  The  results  of  this  experiment  indicate  that  the  basic 
theory  is  valid:  and  that  it  is  complete  in  the  linear  regime. 

This  experiment  underlines  th;  importance  of  extending  the  theory 
to  cover  the  new  generation  of  devices  which  will  operate  in  the 
saturated  or  the  high  sain  regime. 
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p  '  GURE  captions 


Fiauro  1 

The  laser  beam  was  measured  with  a  Reticon  diode  array  to 
determine  its  mode  characteristics.  The  beam  parameter  W(z)  measured  in 
Series  II  is  plotted  and  compared  to  a  TEMoo  Gaussian  fit.  The  beam  is 
well  Focussed  into  the  Klystron*  whose  position  is  indicated  at 
the  center  of  the  horizontal  axis. 


Fisure  2 

Photographs  of  the  transverse  laser  profiles  during  series 
II*  taken  by  the  Reticon  at  two  lonsitudinal  positions!  at  the 
center  of  the  optical  Klystron,  and  4.3  m  in  front  of  it.  The 
horizontal  and  the  vertical  rms  beam  sizes  are  noted.  In  the 
center  of  the  laser  interaction  resion.-  the  laser  profiles  can  be 
fit  well  with  a  Gaussian  shape.  However*  far  from  the  focus*  the 
profiles  acquire  a  significant  non-Gaussian  pedestal  which  is 
arrarent  in  the  photo  on  the  left  (the  discontinuity  in  the 
vertical  Profile  is  an  artifact  caused  by  a  bad  diode). 


Figure  3 


The  electron  beam  measurement  scheme.  Synchrotron  lisht 


omitted  in  the  bending  magnet  before  the  Klystron  is  imaged  onto 
the  Pet  icon. 

Figure  4 

Measured  ratio  of  stimulated  power  to  incident  laser  power 
transmitted  through  an  analyzing  aperture  as  a  function  of  its 
diameter.  The  diameter  of  the  aperture  is  normalized  to  the  laser 
beam  waist  parameter  at  the  measurement  point?  and  the  gain  ratio 
is  normalized  to  unity  at  large  aperture  opening.  The  points  were 
measured  on  two  successive  scans?  the  error  bars  can  be  estimated  by 
the  deviation  of  the  points  from  each  other.  The  solid  curve  is 
the  theoretical  C13  result-  using  the  measured  beam  parameters!  Wo 
=  670  microns  for  the  laser.-  and  <T*~  220  microns  for  the  electron 
beam.  The  dashed  lines  show  the  contributions  of  the  individual 
off-diagonal  terms  in  the  gain  matri:'.-  assuming  that  for  the 
purposes  of  demonstration?  all  the  other  off-diagonal  terms 
are  zero. 

Figure  5 

The  gain  is  measured  as  a  function  of  the  iris  diameter  as  in 
figure  4.-  but  now  for  two  different  electron  beam  sizes.  The 
triangles  taKen  for  the  small  electron  beam  size?  show  the 


largest  deviation  from  unity*  and  therefore  the  largest  excitation 
of  the  hishor  order  modes.  As  explained  in  the  text*  smaller 
electron  beam  sizes  result  in  larser  values  for  the  off-diasonal 
terms  in  the  transverse  sam  matrix.  The  shift  in  the  resonance 
parameter  between  the  two  curves  is  inconsequential:  as  shown  in 
fiaure  G. 

Fiaurc  S 

The  sain  is  measured  as  a  function  of  the  iris  diameter  for  three 
different  values  of  the  dotunins  parameter  which  correspond  to 
a  lin  and  absorption  peaks  displaced  by  four  optical  Klystron 
frinses  (more  than  half  of  the  sain  curve)  from  the  center.  The 
data  all  fall  along  the  same  curve  within  the  error  bars :  and  we 
conclude  that  the  mode  mi  Kina  effect  does  not  depend  on  detun  ms. 
The  changes  in  the  beam  size  are  small*  and  unavoidable  due  to  the 
current  decay  in  ACO  durins  these  measurements. 

Fiaure  7 

The  sain  is  measured  as  a  function  of  aperture  for  a  position 
close  to  the  output  of  the  optical  Klystron.  T h o  error  bars  are 
lar3»  here  because  of  the  constant  real isnmont  of  the  aperture 


position.  The  beam  sizes  are  listed  here  and  in  figure  1.  The 


inr-ut  laser  is  operating  multimode  (see  figure  2).  The  solid  line 
is  the  theorr  for  a  single  mode  input  beam"  me  attribute  the 
dciation  to  the  effects  of  multiple  modes  in  the  input  laser. 


Figure  8 

Two  scans  of  the  gain  us.  diameter  taKen  far  from  the 
interaction  region  under  similar  conditions  to  those  of  figure  7. 
As  was  seen  in  figure  5:  the  scan  with  the  smaller  electron  beam 
(the  solid  points)  shows  a  larger  deviation  from  unit/.  This 
figure  is  different  from  figure  7  because  the  higher  order  modes 
e%  O’  ted  in  the  mode  mixing  experience  a  mode-denendant  phase  shift 
(see  the  text  >  equation  (6)1  in  propagating  awax  from  the 
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